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Abstract: We study matrix quantum mechanics at finite temperature by Monte Carlo 
simulation. The model is obtained by dimensionally reducing lOd U(iV) pure Yang-Mills 
theory to Id. Following Aharony et al., one can view the same model as describing the 
high temperature regime of (l+l)d U(iV) super Yang-Mills theory on a circle. In this inter- 
pretation an analog of the deconfinement transition was conjectured to be a continuation 
of the black- hole/black-string transition in the dual gravity theory. Our detailed analysis 
in the critical regime up to W = 32 suggests the existence of the non-uniform phase, in 
which the eigenvalue distribution of the holonomy matrix is non-uniform but gapless. The 
transition to the gapped phase is of second order. The internal energy is constant (giving 
the ground state energy) in the uniform phase, and rises quadratically in the non-uniform 
phase, which implies that the transition between these two phases is of third order. 
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1. Introduction 

Recently large- N gauge theories are playing increasingly important roles in string theory. 
One of the crucial discoveries was that U(iV) gauge theory appears as a low energy effective 
theory [1] for a stack of N D-branes [2] in string theory. This led to various interesting 
conjectures concerning large-iV gauge theories. For instance, it is conjectured that large- 
N gauge theories obtained by dimensionally reducing lOd U(iV) super Yang-Mills theory 
to D = 0,1,2 dimensions provide non-perturbative formulations of superstring/M theo- 
ries. These are called Matrix theory (D = 1) [3], the IIB matrix model (D = 0) [4] and 
matrix string theory (D = 2) [5], respectively. Another type of conjectures asserts the 
duality between strongly coupled large-iV gauge theory and weakly coupled supergravity. 
In the AdS/CFT correspondence [6-8], for instance, it is conjectured that 4-dimensional 
U(iV) M = 4 super Yang-Mills theory is dual to the type IIB supergravity on AdSs x S 5 . 
Generalizing this correspondence to the finite temperature case, it was argued that the 
deconfinement phase transition on the gauge theory side corresponds to the Hawking-Page 
transition on the gravity side [9]. This has been further extended to non-conformal gauge 
theories including super symmetric matrix quantum mechanics [10]. See refs. [11-20] for 
extensive studies on the relationship between large- N gauge theories at finite temperature 
and the black-hole physics. 

In this paper we study matrix quantum mechanics at finite temperature by Monte 
Carlo simulation. The model is obtained formally by dimensionally reducing lOd U(iV) 
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pure Yang-Mills theory to Id, and it may be viewed as the bosonic part of the low en- 
ergy effective theory of D-particles in type IIA superstring theory [1]. By considering the 
Euclidean time direction as a spatial direction instead, one can view the bosonic model as 
describing the high temperature regime of (l+l)d XJ(N) N = 8 super Yang-Mills theory 
on a circle. In this interpretation an analog of the deconfinement transition was specu- 
lated [15] to be a continuation of the black- hole/black-string phase transition in the dual 
gravity theory 2 . 

Unlike previous works, our detailed analysis in the critical regime up to N = 32 suggests 
the existence of the non-uniform phase, in which the eigenvalue distribution of the holonomy 
matrix is non-uniform but gapless. The transition to the gapped phase appears to be of 
second order. At low temperature, the internal energy is constant (giving the ground state 
energy) as a result of the Eguchi-Kawai equivalence. We use this property to identify 
the uniform phase. As one enters the non-uniform phase increasing the temperature, the 
internal energy starts to rise quadratically. This implies that the transition between the 
uniform phase and the non-uniform phase is of third order. These observations select one 
of the two scenarios suggested by the Landau-Ginzburg analysis [15]. 

The rest of this paper is organized as follows. In section 2 we define the model and the 
observables we study. In section 3 we present an overview of the phase diagram as seen 
from explicit results for the observables. In section 4 we focus on the critical regime and 
study the phase structure in more detail. In section 5 we use our results to speculate on 
the phase diagram of (l+l)d super Yang-Mills theory on a circle. Section 6 is devoted to 
a summary and discussions. In appendix A we derive a formula for the internal energy, 
which is used for its evaluation. In appendix B we present the details of our simulation. 

2. The model and observables 

The model we study in this paper is defined by the action 



where Dt represents the covariant derivative Dt = dt — i [A(t), ■ ]. The one-dimensional 
fields A(t) and Xj(t) (i = 1, 2, ■ ■ ■ ,9) are N x N Hermitian matrices, which can be regarded 
as the gauge field and nine adjoint scalars, respectively. The model can be formally obtained 
by dimensionally reducing lOd U(iV) pure Yang-Mills theory to 1 dimension. The Euclidean 
time t has a finite extent (3, which corresponds to the inverse temperature (3 = 1/T, and 
all the fields obey periodic boundary conditions. 

1 One may also consider the model as the "bosonic Matrix theory" [21], but in that case, the compactified 
Euclidean time direction actually corresponds to a light-cone coordinate, which might make the meaning 
of "finite temperature" a bit subtle unlike in the other interpretations. See refs. [22], however. 

2 Some of the results in refs. [15] were already anticipated in a pioneering work [14] on the phase structure 
of black holes on a circle. There, the relationship of the supergravity solutions to non- and near-extremal 
branes on a circle and to the corresponding dual non-gravitational theories are discussed. 




(2.1) 
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The 't Hooft coupling constant A = g 2 N has the dimension of (mass) 3 , and we fix it 
when we take the large-iV limit. 3 The physical properties of the system depend only on 
the dimensionless effective coupling constant given by 4 

Aeff = p • (2.2) 

In what follows we set A = 1 without loss of generality. 

It is known that the bosonic matrix quantum mechanics undergoes a phase transition 
[15,21,23] at some critical temperature, which can be interpreted as the Hagedorn transition 
in string theory [16,24]. This transition is associated with the spontaneous breakdown of 
the U(l) symmetry 

A(t) h+ A(t) + al , (2.3) 

and therefore it is analogous to the deconfinement transition in ordinary gauge theories. 
The order parameter is given by the Polyakov line 

P = ltvU, (2.4) 

U = Vexp J^dtA{t)^ , (2.5) 

where the symbol Vexp represents the path-ordered exponential and the unitary matrix 
U is called the holonomy matrix. In section 4 we present numerical results suggesting that 
the "deconfined phase" is further divided into two phases. 

As a fundamental quantity in thermodynamics, the free energy T = —^lnZ((3) is 
defined in terms of the partition function 

Z(J3) = j[VX]p[VA}pe- s ^ , (2.6) 

where the suffix of the measure represents the period of the field to be path- integrated. 
However, the free energy T cannot be calculated straightforwardly by Monte Carlo simu- 
lation because that would require evaluation of the partition function Z((3). We therefore 
study the internal energy defined by 

*-£cw (2 ' 7) 

which has equivalent information as the free energy (given that T = E at T = 0). Note 
also that the internal energy at T = provides the ground state energy of the quantum 
mechanical system. In appendix A we derive a formula 

w* E = \ {p2) ' (2 - 8) 

F2 = ~J*L dttT {l X M,Xj(t)?) . (2-9) 



3 This is different from the decompactifying limit [3] in the Matrix theory. 

4 One can confirm this statement by rescaling the fields and the coordinate t appropriately so that all 



the A and T dependence appears in the combination of eq. (2.2). 
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where the symbol ( • ) represents the expectation value with respect to Z{(3). This formula 
enables us to calculate the internal energy E directly by Monte Carlo simulation. As yet 
another quantity, we study 

R2 ^jwL dttr{Xt)2 ' (2,io) 

which represents the "extent of space" . 

The details of our simulation are given in appendix B. The number of sites Nt in the 
Euclidean time direction is taken to be the smallest integer that satisfies N t > Q ^ 2T for 
figs. 1 and 2 , and N t > Q q 5T for figs. 3, 4 and 5. This corresponds to taking the lattice 
spacing to be a ~ 0.02 and a ~ 0.05, respectively. See ref. [25] for a systematic study of 
the finite lattice-spacing effects in an analogous model. Note that the previous work [15] 
used Nt = 5 at any temperature. 

3. An overview of the phase diagram 

In this section we present an overview 
of the phase diagram of the model (2.1). 
Fig. 1 shows the results for the Polyakov 
line. It changes drastically at T ~ 0.9. 
For T > 0.9, the results for N = 16 and 
N = 32 lie on top of each other, and ap- 
proaches the maximum value 1 as T in- 
creases. The results at T > 2 are nicely 
reproduced by the high temperature ex- 
pansion 5 including the next-leading or- 
der terms [26]. For T < 0.9, the ex- 
pectation value (\P\) takes small values, 
which seem to vanish in the large- N limit 
as l/N. 

In fact, the expectation value (|P|) 
in the T — > limit seems to agree well 
with the finite- ./V results obtained by gen- 
erating the holonomy matrix U randomly using the Haar measure, which are represented 
by the horizontal dash-dotted lines in fig. 1 for N = 16, 32. The latter results are seen to 
be proportional to l/N for N = 8, 16, 32, 64. In this way, we can understand the observed 
l/N behavior of (|P|) in the original model at low temperature. On the other hand, finite 
N effects at T > 0.9 seems to be much smaller. This can be understood from the fact that 
in this regime the phenomenological description of the holonomy matrix is given by the 
Gross- Witten model (4.3), which has finite N effects of the order of 0(1/N 2 ) for N S> \/k. 

In fig. 2 we plot the internal energy obtained from the formula (2.8), and the "extent 

of space" (R 2 ). For T < 0.9, we observe that the results are independent of T. This can 

be understood as a consequence of the Eguchi-Kawai equivalence. In general it states that 

the expectation value of a single-trace operator in D-dimensional U(iV) gauge theory is 

5 The results of the high temperature expansion [26] are plotted with the double lines showing the 
statistical errors coming from the Monte Carlo integration over the zero modes. In some cases, the errors 
are so small that one cannot recognize the lines as two separate ones. 




T 



Figure 1: The Polyakov line (\P\) is plotted against 
T for N = 16, and in the critical regime, also for 
N = 32. The dashed line and the solid line represent 
the results of the high temperature expansion [26] for 
N = 16, which are obtained at the leading order and 
at the next-leading order, respectively. 
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Figure 2: The same as fig. 1 but for the internal energy j^E (left) and for the "extent of space" 
(R 2 ) (right). 



independent of the volume in the large- N limit, provided that the U(l) symmetry is not 
spontaneously broken 6 [27]. The T-independence of liniAr^oo -^E at low temperature can 
also be viewed as a consequence of the fact that only U(iV) invariant states appear in the 
low energy spectrum in a confining theory. 



4. Closer look at the critical regime 
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In this section we focus on the critical 
regime and study the phase structure in 
more detail. 

In fig. 3 we plot the Polyakov line 
against temperature magnifying the crit- 
ical regime. For N = 16 the Polyakov 
line changes smoothly with T. However, 
for N = 32 the behavior of the Polyakov 
line changes drastically at T ~ T& = 

0. 905(2). In fig. 4 we plot the internal 
energy E and the "extent of space" (R 2 ) 
in the critical regime. The qualitative 
behavior is similar to the Polyakov line. 

The N = 32 data suggest that all the observables depend continuously on T, but their first 
derivatives with respect to T seem to be discontinuous at T ~ T c \. Thus we conclude that 
there exists a second order phase transition at T ~ T c \ . 

In order to clarify the nature of the phase transition seen above, let us consider the 
eigenvalues of the holonomy matrix U given by eq. (2.5), which we denote as e 1 ® 3 (j = 

1, ■ ■ ■ , N), where ■dj G (—it, it]. Then we define the eigenvalue distribution by 



Figure 3: The Polyakov line (\P\) is plotted against 
T in the critical regime for N = 16 and N = 32. The 
dotted lines are drawn to guide the eye. 



1 N 

7 = 1 



(4.1) 



6 In 3d and 4d pure SU(iV) gauge theory on a torus, the SSB occurs at finite volume [28,29], and the 
relation to the deconfinement transition at finite temperature [30] is discussed. See also refs. [31-33] for 
related works. 
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Figure 4: The same as fig. 3 but for the internal energy E (left) and for the "extent of space" (R 2 ) 
(right). The solid line in the left panel is a fit to — eo = c(T — T C 2) P , where eo is given by the 
low-temperature data. An analogous fit in the right panel yields consistent values for T C 2 and p. 
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Figure 5: (Left) The eigenvalue distribution p{6) is plotted for N = 32 at three different values 
of T in the critical regime. The solid lines represent a fit to the Gross- Witten form given by either 
(4.4) or (4.5) depending on the fitting parameter k. (Right) The value of k obtained by fitting the 
eigenvalue distribution to the Gross- Witten form is plotted against T. The dotted line is drawn to 
guide the eye. 

where we assume that the U(l) transformation (2.3) is applied to each configuration in 
such a way that tr U becomes real positive before taking the expectation value. 7 Using this 
definition, we have 

(\P\)= d9p(8)e* 6 . (4.2) 

J — IT 

Fig. 5 (left) shows the result of the eigenvalue distribution p{0) for N = 32. It is clear 
that the second order phase transition is associated with the emergence of a gap in the 
eigenvalue distribution. 

7 Note that the U(l) transformation (2.3) rotates all the eigenvalues e^ by some constant angle on 
the complex plane. If we took the expectation value in eq. (4.1) naively, we would trivially get a uniform 
distribution. 
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Thus we find that the "deconfined phase", in which the U(l) symmetry (2.3) is spon- 
taneously broken, is further divided into the gapped phase and the non-uniform phase. In 
these phases, we may consider the Gross- Witten model [34] 

Z GW = J dUexp^(trU + trU^)\ (4.3) 

as a phenomenological model for the holonomy matrix, where the parameter k should be 
determined as a function of T. The large-iV limit of the eigenvalue distribution for the 
Gross- Witten model, which we denote as Pgw(6), is known analytically [34]. For k < 2 we 
have a gapped distribution given by 



PGwW = ^(cos^W|-sin^ (4.4) 



for |0| < 2 sin 1 and otherwise. For n > 2 we have a gapless distribution 



PGw(0) = ^(l + ^cos0) ( L.5) 



From these results, one obtains [34] 
1 



Jim / J_ tr u) = r d 6 PGW (6) e» = \\ 1 for K < 2 > (4.6) 
iV-oo\iV / GW J_ n ^ GWV ; I I for K >2, 

which crosses 1/2 at the critical point n = 2. Note that (4.6) and its first derivative with 
respect to k is continuous at k = 2, but the second derivative has a discontinuity. Thus, 
the Gross- Witten model undergoes a third order phase transition at n = 2. 

As we see from fig. 5 (left), the distribution p{9) for N = 32 can be nicely fitted to 
the Gross- Witten form by choosing k appropriately at each T. The value of k obtained in 
this way is plotted against T on the right panel of the same figure. We observe that the 
first derivative of k with respect to T is discontinuous at T = T c ±. This is reflected in the 
behavior of the Polyakov line (\P\) for N = 32 shown in fig. 3. Thus it is possible to obtain 
a second (instead of third) order phase transition between the non-uniform phase and the 
gapped phase in the present model, despite the fact that the eigenvalue distribution is well 
described by the Gross- Witten form. 

In the "confined phase", the U(l) symmetry (2.3) is unbroken, and therefore the eigen- 
value distribution p{9) is uniform in the large- iV limit; hence we call it the uniform phase 
following the present terminology. However, as we have seen in fig. 1, the Polyakov line, 
which should vanish for a uniform distribution, seems to be of 0(1/N) at low temperature, 
which actually vanishes slowly with increasing N. On the other hand, the Eguchi-Kawai 
equivalence, which is a consequence of the unbroken U(l) symmetry, holds with high accu- 
racy at low temperature as we have seen in fig. 2. In general the breaking of the Eguchi- 
Kawai equivalence due to finite N effects is expected to be of 0(1/A^ 2 ). Taking advantage 
of this fact, we identify the uniform phase using the Eguchi-Kawai equivalence instead of 
identifying it directly using the eigenvalue distribution. In particular, the internal energy 



-7- 



normalized by iV 2 is constant at low temperature, and the constant value, which gives the 
ground state energy, is extracted as £o = 6.695(5) from the N = 32 data. With increasing 
T, one enters the non-uniform phase at some critical point T C 2, where the internal energy 
starts to deviate from this constant value. We fit the results of the internal energy for 
N = 32 shown in fig. 4 to the behavior 



with three free parameters, which are determined as c = 413 ± 310, T C 2 = 0.8758(9) and 
p = 2.1(2). We redo similar analysis for the "extent of space" (R 2 ). The constant value at 
low temperature is given by (r ) 2 = 2.291(1). The deviation (i? 2 )-(r ) 2 can be fitted to the 
behavior (4.7) with three free parameters, which are determined this time as c = 39 ± 36, 
T C 2 = 0.8763(4) andp = 1.9(2). The results for T C 2 and p obtained from the two observables 
are consistent with each other. Moreover they suggest p = 2, which implies that the phase 
transition between the uniform phase and the non-uniform phase is of third order. 

Thus we have identified the non-uniform phase between the two critical points T c \ = 
0.905(2) and T C 2 = 0.8761(3). The range of T is very narrow. This is not so surprising, 
however, given that the Polyakov line grows very rapidly as we increase the temperature 
in the critical regime. Note that the Polyakov line should lie within the range [0, \] m the 
non-uniform phase according to the (successful) phenomenological description in terms of 
the Gross- Witten model. 

Let us also comment that the phase structure obtained above is consistent with one of 
the two scenarios suggested by the Landau-Ginzburg analysis [15]. The other scenario was 
that the Polyakov line jumps from to 1/2 at some critical temperature, indicating a single 
first order transition between the uniform phase and the gapped phase. This behavior was 
observed in the plane- wave matrix model at finite temperature [35,36]. On the other hand, 
in ref. [23] we observed no phase transition in the fuzzy sphere background. This might 
also be the case with the super symmetric version of the present model [11, 18]. 

5. Connection to the black-hole/black-string transition 

In fact one can show that the model (2.1) describes the high temperature regime of (l+l)d 
super Yang-Mills theory on a circle. At low-temperature strong-coupling regime, the 2d 
theory is expected to have a first order phase transition, which corresponds to the black- 
hole/black-string transition in the dual gravity description. In this section we review this 
connection, and discuss the implication of our results on it. 

Let us consider 2d U(iV) J\f = 8 super Yang-Mills theory, 8 which can be obtained by 
dimensionally reducing lOd U(iV) M = 1 super Yang-Mills theory to 2d. In order to put the 
theory at finite temperature, we compactify the Euclidean time direction to a circle with 
the circumference (3, which corresponds to the inverse temperature (3 = i. Furthermore 

8 We put tildes on all the parameters of the 2d theory to distinguish them from the parameters of our 
Id theory. 
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Figure 6: A schematic view of the phase diagram of the 2d super Yang- Mills theory in terms 
of dimcnsionless parameters TL and XL 2 . The region below the dotted line (5.1) can be well 
described by our Id model, from which we obtain the critical lines corresponding to (5.4). The 
upper left corner is conjectured to have a dual gravity description, which predicts the first order 
phase transition at (5.5). The dashed lines represent our speculation that the non-uniform phase 
ceases to exist at a tri-critical point for consistency with the prediction from the gauge/gravity 
duality. 

we compactify the spatial direction to a circle 9 with the circumference L. At sufficiently 
high temperature, the temporal Kaluza-Klein modes decouple, and one obtains our model 
(2.1) with appropriate identification of parameters. Note in particular that the fermions in 
the 2d theory decouple due to anti-periodic boundary conditions in the temporal direction, 
and therefore one ends up with a bosonic theory. The condition for the decoupling of the 
temporal Kaluza-Klein modes is given by 

f>> {ij' 3 - <51) 

In this interpretation the i-direction of our Id model is identified with the spatial 
direction of the original 2d theory, and the relationship among the parameters is given by 

p = L, A = AT . (5.2) 

Hence the dimensionless effective coupling constant X e g defined in eq. (2.2) is written as 

Aefr = ATX 3 . (5.3) 



9 This model is formally the same as the matrix string theory [5,37,38] except that fermions obey anti- 
periodic boundary conditions in the temporal direction. 
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Note also that the Polyakov line (2.4) should be regarded as the Wilson loop winding 
around a spatial direction in the 2d theory. In ref. [15] the Id model (2.1) was studied 
from this point of view and a phase transition was observed around \ c g ~ 1.4. Our results 
suggest that actually there exist two phase transitions at 

A cfr ~ 1.35(1) and A cff ~ 1.487(2) , (5.4) 

which are of second order and of third order, respectively. 

At large A and small T, the 2d theory has a dual gravity description, which predicts a 
first order phase transition at [15] 

TL = . (5.5) 

VAL 2 

On the gravity side, this corresponds to the black-hole/black-string transition, which is 
associated with the Gregory-Laflamme instability [39] of the black string winding around 
the spatial direction. The black-string phase and the black-hole phase can be naturally 
identified with the uniform phase and the gapped phase on the gauge theory side. A similar 
correspondence was suggested earlier by ref. [40]. 

We speculate that the first order phase transition predicted at low temperature actually 
splits into two continuous transitions as one increases the temperature. In figure 6 we 
present a schematic view of the phase diagram of the 2d super Yang-Mills theory that 
emerges from the present work. 



6. Summary and discussions 

In this paper we have investigated the phase structure of matrix quantum mechanics at 
finite temperature. We have identified three phases. At high temperature, the high tem- 
perature expansion up to the next-leading order gives a precise description of various ob- 
servables. At low temperature, the internal energy, in particular, does not depend on the 
temperature as a consequence of the Eguchi-Kawai equivalence. This property enables us 
to determine the critical point and the order of the transition between the uniform phase 
and the non-uniform phase. 

In the non-uniform phase and the gapped phase, the eigenvalue distribution of the 
holonomy matrix can be fitted nicely to the Gross- Witten form. While this suggests that 
we are already seeing the large- N behaviors, it is certainly desirable to confirm the stability 
of our results against increasing N, which we leave for future investigations. 

Our results can be alternatively interpreted as representing the high temperature be- 
havior of 2d N = 8 super Yang-Mills theory. The low temperature behavior of that theory 
can be predicted by the gauge/gravity correspondence. We speculate that the non-uniform 
phase identified in the present paper ceases to exist below some temperature in the phase 
diagram of the 2d super Yang-Mills theory. If so, it is interesting to locate this point 
explicitly. 

The non-uniform phase may exist also in higher dimensional bosonic gauge theories on 
a finite torus [28,29]. That will have implications on the phase diagram of super Yang-Mills 
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theories in D = 3, 4. The low temperature regime of these theories is discussed in ref. [42] 
based on the gauge/gravity correspondence. 

It would be also interesting to study the supersymmetric version of the present model 
using the non-lattice simulation proposed recently [43]. In particular, by studying the 
low temperature regime, we would be able to test the predictions of the gauge/gravity 
correspondence directly. The phase transitions are expected to disappear [11,18], and the 
internal energy is expected to vanish as T — > with the power law behavior [13] obtained 
from the dual black-hole geometry [41]. Indeed the preliminary results for a simplified 
model with 4 supercharges [43] agree qualitatively with these expectations. Studies of the 
model with 16 supercharges are in progress [44]. 
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A. Derivation of a formula for the internal energy 

In this section we derive the formula (2.8) relating the internal energy of the present model 
to the expectation value (2.9), which is directly accessible by Monte Carlo simulation. 
Let us first rewrite (2.7) as 

E = h fbl, (A.1) 

Z(J3) a/3^o A/3 y ' 

where (3' = (3 + A/3, and represent Z{(3') for later convenience as 

Z((3') = J[VX'] p ,[VA%,e- s ' , (A.2) 

where S' is obtained from S given in (2.1) by replacing /3, t, A(t), Xi{t) with f3', t', A'(t'), 
X-(t'). In order to relate Z([3') to Z((3), we consider the transformation 

t' = ^t, A'(t') = lA(t), X' t (t') = ^ X t (t) . (A.3) 

The factors in front of the fields are motivated on dimensional grounds, and in particular 
we have \DX'\^ = \DX\p and [DA']gi = [DA]/}. Under this transformation, the kinetic 
term in S' reduces to that in S, but the interaction term transforms non-trivially as 

jfVtr ([X'it'lXW)} 2 ) = (|) 3 j^dttr ([*(*), A>(f)] 2 ) . (A.4) 

This gives us the relation 

Z(/3') = Z(/3)|l-|iV 2 A/3(F 2 ) + 0((A/3) 2 )| , (A.5) 

where the operator F 2 is defined by (2.9). Plugging this into (A.l), we get (2.8). 
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B. Details of Monte Carlo simulation 

In this section we present the details of our Monte Carlo simulation. We discretize the 
Euclidean time direction and obtain the partition function 

Zlat = yVp*i] exp(-Slat) (B.l) 

4, . „*|>{i C / <">-^" +1 fW-^-W ) 2 - i W („).^(„)p} , (B.2, 

where a is the lattice spacing, and the inverse temperature is given by (3 = aN t . The link 
variables V{n) are N x N unitary matrices representing the gauge connection between the 
lattice sites. Due to the periodic boundary conditions, we have Xi{Nt + 1) = -Xj(l). 

Although it is possible to simulate the system (B.l) directly, let us simplify it [15] by 
taking the static diagonal gauge 

V(l) = V{2) = ■■■ = V(N t ) = V = di3,g(e iei/Nt ,e W2/Nt ,- ■ ■ ,e^ /7Vt ) . (B.3) 
The integration measure for the angular variables 6 a G (— ir,7r] is given by 

A(0) = J] S m 2 (^) , (B.5) 



where A(0) is the Vandermonde determinant. 

The operators (2.5), (2.9) and (2.10) can be calculated on the lattice by 

U = V(l)V(2)---V(N t ) =diag(V\e* e V-- ,e Wjv ) , (B.6) 
1 Nt 

p2 = -^E tr ( [Xi(n) '^ (n)]2 ) ' (R7) 

* n=l 
1 Nt 

R2 = m E^{ x ^) 2 )- ( B -8) 



1 n=l 



In order to employ the heat-bath algorithm for updating -Xj(n), we use the trick used in 
simulating the bosonic IKKT model [45]. Namely, we introduce N x N hermitian matrices 
Qij(n) = 1, ■ ■ ■ , 9; i < j) as auxiliary fields, and consider the partition function 

Zte = j[de][dX}[dQ}exp(-S Ut ) (B.9) 
1 (VXi(n + l)Vi - Xi(n) s 2 



Aft 

Si at = aN ^2 tr 



n=l 



+^E(^( n ) 2 - 2 ^( n )^( n )>^'( n )}) + 2 E^( n ) 2 ^N 



(B.10) 
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Integrating out the auxiliary fields Qij(n), we retrieve the original action (B.2). We update 
Qij(n) and Xi(n) for each n in the same way as described in ref. [45]. 

After updating all the elements of Qij(n) and Xj(n) for each n, we update the angular 
variables 6 a using the standard Metropolis algorithm. Typically the acceptance rate is 
not very high. We therefore repeat the Metropolis procedure sufficiently many times so 
that most of 6 a get updated. This defines our "one sweep". For N = 16 (N = 32) 
we have made 200,000 sweeps (150,000 sweeps) in total, and discarded the first 20,000 
sweeps (40,000 sweeps) for thermalization at each temperature. The simulation has been 
performed on PCs with Pentium 4 (3GHz), and it took a few weeks to get results at each 
temperature for N = 32. 
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